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Abstract: This paper analyzes the eigenvalue distribution of neutral differential systems and the cor- 
responding difference systems, and establishes the relationship between the eigenvalue distribution and 
delay-independent stability of neutral differential systems. By using the "Complete Discrimination Sys- 
tem for Polynomials" , easily testable necessary and sufficient algebraic criteria for delay-independent 
stability of a class of neutral differential systems are established. The algebraic criteria generalize and 
unify the relevant results in the literature. Moreover, the maximal delay bound guaranteeing stability 
can be determined if the systems are not delay-independent stable. Some numerical examples are 
provided to illustrate the effectiveness of our results. 
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Introduction 



Neutral differential system model can often be found in control process, physics, chemical engineer- 
ing and ecology [3 El 1101 1^ • Stabihty of neutral differential systems is an important performance 
index and has been investigated by many authors. From the control-theoretic viewpoint, the 
delay-independent stability of a delay system ensures robustness and reliabiUty of the system for 
all delays. It is well known that the asymptotic behavior of the zero solution of a linear autonomous 
retarded differential system is determined by its eigenvalues. Just as ordinary differential systems, 
the stability of the zero solution is equivalent to all eigenvalues being with negative real parts 
121 El El El El El El EOl • However, it is more complex when considering the neutral differential 
system, i.e., all eigenvalues being with negative real parts are not sufficient for the stability of 
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the zero solution pilHI51l^ ll(J[[T^I21| . and the global hyperbolicity of a corresponding difference 
system must be taken into account ^ ^ . In neutral differential system, when characterizing 
system stability by means of the corresponding characteristic equation, it is usually required that 
the supremum of the real parts of all eigenvalues is negative 13 13 UHl CSl EH , which is notably 
different from retarded differential system. 

Consider the linear autonomous neutral differential system 

N N 

x{t) - J2 Bkiit - 7fc • r) = Aox{t) + ^ Akx{t - jk ■ r) (1) 
fc=i fe=i 

where i?+ = [0,oo),i?= (-oo, +oo), x G i?", G i?"^", e G i?"X",r = (n,- •■,rM) G 

(i?+)^^7fc = (7fci,---,7fcM).7fc, > are integers, 7^ ^ 0,7^ -r = J2jLilk,rj,k = l,2,---,iV,j = 
1, 2, • • • , M. Suppose det[/ — J2k=i ^ (the notation det[-] stands for the determinant of the 
corresponding matrix). The characteristic equation of System (1) is 

N N 

G(A, r. A, B) = det[A(/ - ^ Bke-^'"--'') - - ^ Afee"^'^'^ ''] = (2) 

k=l k=l 

where / stands for the n x n identity matrix, A = {Aq, Ai, • ■ • , An), B = {Bi, • • • , -Bat). 

Definition 1 [H] System is said to be delay-independent stable, if Vr G (i?^)^^, the 
supremum of the real part of A satisfying Equation ^ is negative, that is, there exists a (5 > , 
such that {ReA : G(A, r. A, B) = 0} (1 [-6, 00) = <j). 

In the literature, delay-independent stability is also called all-delay stability, unconditional 
stability or absolute stability. 

Consider the following difference system which is closely related to System |Q 

N 

x{t) - ^ Bkx{t - 7fc • r) = (3) 
fe=i 

The characteristic equation of System Q is 

N 

£;(A,r,S) =det[/-^Bfee-^^'=n = (4) 

k=l 

Definition 2 [5] System © is said to be globally hyperbolic at B, if Vr G the 

real part of A satisfying Equation Q does not intersect with a neighborhood of zero, that is, there 
exists a 5 > 0, such that {ReA : E{\ r, B) = 0} Ci [-S, S] = if). 

In the papers pi llUlIT^ . some elegant theoretical results on the stability of the neutral differ- 
ential system have been established. One important result is 

Lemma 1 !^ Let 

N N 

p(A,si,---,SM,AS) =det[A/-(/-^i?fcS^ + •••^m")] 

fc=i 

then. System is delay-independent stability if and only if 

(HI) System y{t) — ^^=1 Bky{t — 7fc • r) = is globally hyperbolic at B] 

{H2) ReA[(/ — X^fcLi Bk)~^ SfcLo ^k] < 0, where ReA[-] denotes the real part of eigenvalues of 
the corresponding matrix; 

(HS) p{iy,Si,---,SM,A,B)^0,yyeR-{0},\sj\ = l,j = l,2,---,M . 
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The above lemma gives necessary and sufHcient conditions on delay-independent stability of a 
general linear autonomous neutral differential system, which is the fundamental basis for studying 
neutral differential systems. However, it is difficult to apply these conditions in practice, because 
{HI) and {H3) are both transcendental and can't be explicitly tested in engineering computations. 
Therefore, many authors in the area of differential systems tried to derive simple easily-testable 
stability criteria from various viewpoints |11II13[[TH1|19[ 1^1221 1231 124[I25| . These papers considered 
some special systems and just obtained some sufficient conditions, furthermore, most of which have 
never considered the global hyperbolicity of the corresponding difference systems. 

In this paper, we consider a more general linear multi-delay neutral differential system as follows 

N N 

x{t) - Bkx{t ~ kr) = Aox{t) + Akx{t - kr) (6) 

fe=l k=l 

where x e i?",Ao G R"''",Ak G R'"''\Bk £ i?"''",T e R+,k = 1,2, •■•,7V, and we assume 
det[/ — ^k] 7^ 0. Based on Lemma 1, and discussions on the global hyperbolicity of the 

corresponding difference system QJE], algebraic criteria for delay- independent stability of the above 
system will be established by means of the "Complete Discrimination System for Polynomials" 
\27i I28L I29j . The algebraic criteria generalize and unify the relevant results in the literature. 
Moreover, the maximal delay bound guaranteeing stability can also be determined when the system 
is not delay-independent stable. 

This paper is arranged as follows: In section 2, the relevant results of the "Complete Dis- 
crimination System for Polynomials" are presented; Section 3 discusses the global hyperbolicity 
of difference system y{t) — J2k=i ^kvit — kr) = 0; Section 4 analyzes the delay-independent sta- 
bility of System © and establishes the easily testable algebraic criteria; Section 5 presents some 
corollaries and examples by means of the algebraic criteria provided in section 4; Finally, a simple 
conclusion is given in section 6. 



2 Complete Discrimination System for Polynomials 

From Lemma 1 in section 1, the problem of delay- independent stability of System ((^J), boils down 
to analyzing the eigenvalues' distribution of System © and the corresponding difference system. 
In the sequel, we will transform this problem into the real root distribution of some related poly- 
nomials. The following is some results of the "Complete Discrimination System for Polynomials" , 
which can effectively discriminate the distribution of the roots of polynomials by a series of explicit 
expressions and calculations with polynomial coefficients [HI El UHl HSI • 
Suppose 

/(x) = aox" -l-aix"-i + ha„ (7) 

the Sylvester matrix 0121112111221 of f{x) and its derivative f'{x) is defined as 

ao fli a2 
naQ in — l)ai 
ao ai 
noo 

ao ai • • • an 
nao • ■ • a„_i 



2a„_2 fln-i 
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which is caUed the discrimination matrix of f{x), denoted as Discr{f). 



[Di{f),D2{f ),■■■, DM)] 

is a sequence of the determinants of the first n even-order principal sub-matrixes of Discr{f), 
formed by the first 2k rows and first 2k columns, k — 1,2, n, which is called the discrimination 
sequence of f{x). Furthermore 

[sign{Di),sign{D2), ■ ■■ , sign(Dn)] 

is called the sign list of the discrimination sequence[Z?i, • • • , £)„], where sign{-) is the sign 
function, namely 



1 


if X 


>o, 





if X 


-0, 


-1 


if X 


< 0. 



Construct a revised sign list [ei,e2, - ■ ■ , for a given sign list [si, S2, • • ■ , s„] as follows: 
1) If [si, Si+i, ■ • ■ , Si+j] is a section of the given list, where Si ^ 0; Si+i = = ■ ■ • = Si^j^i 
0; Sij^j 7^ 0, then replace the subsection 



[Si+l, Si+2, • • • , Si+j-i] 



by 



r -H, 



namely, Si+r = (^1) ^ ■ Si,r = 1,2, - ■■, j — \, where [a] denotes the greatest integer equal to 
or smaller than a. 

2) Otherwise, let £k = Sk, i-e. no change is made for other terms. 

Lemma 2 [231211 ES| Given a real polynomial f{x) = aox" + aix""^ H h a„ e P", if the 

number of sign changes in the revised sign list of its discrimination sequence is and the number 
of nonzero elements in the revised sign list of its discriminant sequence is fi, then the number of 
the distinct real roots of f{x) is /i — 2j/. 

Remark 1 The discrimination sequence of polynomial f{x) can also be constructed by the 
principal sub-matrices of Bezout matrix [2111211 of f{x) and f'{x); the number of distinct real roots 
of the polynomial /(x) can also be determined by the sign differences of Bezout matrix [2Hl[2n| of 
f{x) and f{x) . 

Remark 2 The original complete discrimination system of polynomial [23 [2H1 [2S1 is more 
general than Lemma 2, which can also be used to determine the number of complex roots and the 
multiplicities of repeated roots. 



3 The Global Hyperbolicity of Difference Systems 

Consider the difference system 



N 

x['] 

fc=l 

and its characteristic equation 



(t) - ^ Bkx{t - 7fc ■ = (8) 
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N 

E{X,r,B) =det[I ~Y,Bke"^^'' ''] = (9) 

k=l 

where all parameters are the same as before. In order to analyze the global hyperbolicity of the 
difference system (jSJ, for its characteristic equation by we have 

Lemma 3 1 Suppose ri, r2, • • • , tm are commensurable, that is, there exist a /3 > and some 
integers rij, such that rj, — k — 1, • • • , M, then there exists an integer p, such that E{X, r, B) 
is a polynomial of some degree p in e~^^ . Denote the p'^'-degree coefhcient of e^"^^ as Ap, then 

p 

E{X,r,B)=ApY[{e-^^ -s,) 

and 

Z{B, r) = Z{B, r) = {-^ In |,s,|, = 1, 2, • • • ,p} 

P 

where Z{B, r) = {ReA : E{\ r, B) = 0}, r) = Z(B,r) = cl{Z(B, r)} is the closure of Z[B, r). 
Now, consider the corresponding difference system of System © 

N 

y{t)-Y,Bky{t-kT)^0 (10) 
fc=i 

Its characteristic equation is 

AT 

e{X,T,B) = det[I-Y,Bke-^''l = (11) 
Lemma 4 The difference system Hl()|l is globally hyperbolic at B if and only if 

N 

\/9 e [0, 27r], det[/ - ^ Bke'''''] ^ 0. 
fe=i 

Proof. From Lemma 3, the number of the roots of the characteristic equation (|11|) is finite, 
and distribute discretely in R x R plane. So the following conclusion comes naturally: 

If the real part of all characteristic roots is nonzero, then there exists a 5 > 0, such that 
{ReA : E{X, r, B) = 0} n [-S, S] = 0. 

That is, the difference system (|1U|I is globally hyperbolic at B if and only if 

N 

Vy e R, e{iy, r, B) = det[I - ^ BkC-'y''^] ^ 0. 

fe=i 

Letting — yr = 6, the lemma is proved. □ 

Owing to the above lemma, it may be feasible to transform the problem of global hyperbolicity 
of the difference system ()1(J|I into the real root existence of some related polynomials. 

4 Algebraic Criteria for Delay- independent Stability of Neu- 
tral Differential Systems 

Consider the linear multi-delay neutral differential system 

N N 

x(t) - ^ Bkx{t ~ kr) = A(^x{t) + ^ Akx{t - fcr) (12) 

k=l k=l 
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The characteristic equation is 

N N 



H{X, r, A, B) = dct[A(/ - ^ Bfce~^'^^) - - ^ ^^6"^'=^ = (13) 



fc=i fc=i 

The corresponding difference system is 

N 



y{t)-J2Bky{t-kT)=0 (14) 



k=l 

and its characteristic equation is 

JV 



e(A,r,B) = det[/-^Sfce-^'=^] = (15) 
Theorem 1 The neutral differential system (|12|l is delay-independent stable if and only if 

(i) det[/ - Ef=i Bke^'^O] ^0,y9e [0, 2^], 

(ii) ReA[(/ - Eti Bk)-' Eto ^'^] < 0, 

(iii) det[iy(/ - Bke'"') ~ Ao ^ ^fee^'l 0, V0 e [0, 27r], Vj; G i?, y ^ 0. 
Proof. By lemma 1 and lemma 4, letting —yr = 6, the theorem is proved. □ 

Condition (i) and condition (iii) are also transcendental, which are still difficult to test numer- 
ically. Next, we will give a simple criterion. 

1 + z 

The bilinear transform u; = is a one-to-one mapping between the set {w = e** : 6* G [0, 2tt]} 

I — z 

and the set {z — iy : y E R}. By this transformation, condition (i) in Theorem 1 is equivalent to 

det[/-^i?.(i±^)'=]^0 
fc=i 



that is 



Let 



AT 



det[(l - izfl - ^ Bfc(l + iz)\l - izf-''] ^ 



fc=i 



N 
k=l 



d{z) = det[(l - iz)^I - Bk{l + izf{l 

k=l 

and f{z) — Re[(i(z)], (7(z) — lm[d{z)], where f{z),g{z) are real polynomials in z. Suppose 

I /(z) = aoz'i -l-aiz'i-i H hazj, 

\ g{z) = hoz^' + + • • • + 

where ^ 0, 6o 0. 

Then condition (i) in Theorem 1 is equivalent to the condition that pt)|l has no common real 
roots. 

By the same transformation, condition (iii) in Theorem 1 is equivalent to 

N N 

^ — ' 1 ~ iz ^ — ' 1 — iz 

k=l k=l 
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namely, 

det[(zy/ - ^o)(l - izf - + iz)\l - iz)^-'' - ^ ^,.(1 + iz)\l - iz)""-^] ^ 

Let 

D{z, y) = det[(zy/ - A„){1 - iz)"" - iyY,Bk{l + izf(l - iz)^-'= - ^ Ak{l + izf(l - iz) 



AT 



fc=l 



N 



N 



k=l 



N 



N-ki 



fc=l 



k=l 



(17) 



and F{z,y) = Iic[D{z,y)],G{z,y) = lni[D{z,y)], where F{z,y),G{z,y) are real polynomials in 
{z,y). Suppose 

f Fiz, y) = ao{y)z^^ + ai(y)z'^-i + • • • + a;, (y), 

\ G(z,2;) = fooCy)^"^ + h^{y)z^--^ + ■■■ + h^M- 

Then condition (iii) in Theorem 1 is equivalent to the condition that p7l) has no common real 

roots (z, y),z e R,y e R~ {0}. 

Before presenting the stability criteria, we first introduce two lemmas in algebra theory EHl • 
Lemma 5 The equation (|16|l has at least one common root if and only if its resultant equals 

zero, i.e.. 



= det 



ao ai a2 
ao ai 



6o hi 62 
60 bi 



■■■ ail 



ao ai 

■ ■ ■ hm. 



ail 



ail 



> mi rows 



1 1 rows 



= 0. 



(18) 



Lemma 6 Suppose ao(?/) or 6o(y) 7^ 0, the equation H17() has at least one common root 
(z, y) and y G i? if and only if its resultant 



R{F, G) = det 



ao(y) ai(y) a2(y) 
ao(y) ai(y) 



ahiv) 



ahiy) 



ao{y) ai(y) 

&o(y) bi{y) &2(y) •■• &m2(y) 

^o(y) bi{y) ■■■ brnAv) 



My) bi{y) 



ahiy) 



bm2iy) 



> 1712 rows 



I2 rows 



has a real root y E R. 

Similarly, F{z,y), G{z,y) can be rewritten in the following formula: 

f F(z, y) = co(z)y'3 + ci{z)y^^-^ + ■■■+ ci,{z), 
I G(z, y) = do(z)y™^ + di(z)y™^-i + ■ • • + 4„3(z). 



(19) 



(20) 



Then condition (iii) in Theorem 1 is equivalent to the condition that ()2U|I has no common real 
roots (z, y), z G -R, y G -R — {0}. 

Also from algebra theory |H1 EH] , we have: 
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Lemma 7 Suppose co{z) / or do{z) ^ 0, the equation (|2(J|I has at least one common root 
(z, y) and z £ i? if and only if its resultant 



R{F, G) = det 



co{z) Ci(z) C2{z) 

co{z) ci{z) 



do{z) di{z) d2{z) 
dQ{z) di{z) 



ci^{z) 



co{z) ci{z) 

' ' ' dm^ (^) 



ci.iz) 



do{z) di{z) ■■■ d„^^{z) 



ma rows 



> h rows 



(21) 



has a real root z & R. 

From Lemma 6 and Lemma 7, we have: 

Lemma 8 If {zo,yo) are a pair of real roots of F{z,y) — 0,G{z,y) — 0, then R{F,G){yo) = 
0,i?(F,G)(zo) = 0. 

Proof. The lemma is proved by using lemma 6 and lemma 7 directly. □ 

Remark 4 It should be noticed that the converse proposition of Lemma 8 is not correct, which 
can be tested by some examples easily. 

Now we are in position to present our main result. The algebraic criteria for the delay- 
independent stability of the neutral differential system (|12ll are as follows: 

Theorem 2 The neutral differential system (|12|l is delay-independent stable if and only if 

(i) 116(1 has no common real roots. 

(ii) RcA[(/ - Eti Bk)-' Eto M < 0. 

(iii) H17|l for t^'20\i ) has no common real roots (z, y) , z ^ R, y ^ R — {0}. 

Proof. From the above discussion, condition (i) in theorem 1 is equivalent to ((16|l having no 
common real roots; condition (ii) in theorem 2 is same as condition (ii) in theorem 1; condition (iii) 
in theorem 1 is equivalent to H17|l (or (|20|) 'l having no common real roots {z , y) , z G R, y G R — {0}. 
This completes the proof. □ 

By Lemma 5, the condition (i) above is equivalent to one of the following conditions: 

(al) the determinant R{f,g) ^ 0. 

(a2) R{f,g) — and g{z) ^ 0, where z is the real root of /(z), and R{f,g) is the resultant of 
(Unj). 

(a3) R{f,g) — and f{z) ^ 0, where z is the real root of g{z), and R{f,g) is the resultant of 

m- 

Similarly, by Lemma 6, the condition (iii) above is equivalent to one of the following conditions: 
(bl) when ao{y) 7^ or bo{y) ^ 0, R{f,g) has no nonzero real root, where R{f,g) is the 
resultant of lfT7|) . 

(b2) when ao{y) 7^ or 60(2/) 7^ 0, R{f,g) has nonzero real root. But when taking the root 
back into ((17|l . then H17() has no common real root. 

(b3) when there exists a real number y, such that ao(j/) = and &o(j/) = 0, then take this real 
number y into H17|l . and H17|l has no common real root. 

Dually, by Lemma 7, the condition (iii) above is also equivalent to one of the following condi- 
tions: 

(cl) when co{z) ^ or do{z) ^ 0, -R(/, (?) has no nonzero real root, where i?(/, g) is the resultant 
of (EOJ. 
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(c2) when co{z) / or do{z) ^ 0, -R(/, g) has nonzero real root. But when taking the root back 
into (|20|l . then H2U|I has no common real root. 

(c3) when there exists a real number z, such that c^i^z) — and d^i^z) = 0, then take this real 
number z into (|20|l . and (|20|l has no common real root. 

Furthermore, by Lemma 8, the condition (iii) above can also be checked with the following 
condition: 

(dl) Vyo e {v\R(F, G){y) = 0} n {i? - {0}}, Vzq G {z\R{F, G){z) = 0} n i?, (yo, zo) are not the 
roots of Equations (|17|l or Equations (|20() . 

Remark 5 The conditions in Theorem 2 are all algebraic conditions. Condition (ii) can 
be checked by the well-known Hurwitz Criterion [H]. Determining real roots of polynomials in 
conditions (i) and (iii) in Theorem 2 can be carried out by the "Complete Discrimination System 
for Polynomials" mentioned before. 

Remark 6 The conditions in Theorem 2 are necessary and sufficient. The sign list of 
the discrimination sequence of polynomials with symbolic coefficients can be obtained easily by 
computer |27l 1281 1^ . Therefore, "on-line" determining the delay-independent stability of neutral 
differential systems l)12f) can be realized, namely, an efficient algorithm can be set up by this 
theorem. 

Remark 7 System H12II is very general, which covers various forms of systems studied in 
ISliainiinillSliniliniEniESESlElESlESllSIl- Theorem l and Theorem 2 generalize the 

relevant results in |71 EO]. More specifically, if Bfc = 0, fc = 1, • • • , iV, in (jHI, then System lO 
degenerates to a retarded differential system. The results in this paper are completely consistent 
with that in 01301 ■ 

Just like the discussion in "7", when System l|12() is not delay-independent stable, the maximal 
delay bound guaranteeing stability can be determined as following 

T:=min|r|r>0,(z,y)e I ^['J^J^Zq , ^ £ i?, Z/ e i? - {0}| , e'^ = i-^^, -yr = . (22) 

In the following section, some numerical examples are provided to to test the delay-independent 
stability of delay systems and compute the maximal delay bound when the systems are not delay- 
independent stable. 



5 Corollaries and Examples 

In what follows, we will present some explicit algebraic criteria for some simple neutral differential 
systems. 

Corollary 1 The system 

x{t) + cx{t - r) + ax{t) + bx{t - r) = (l-f c ^ 0) (23) 

is delay-independent stable if and only if 

(1 - c){b-a) < and {l + c){b + a) > 

Proof. Let 

f{z) = Re[(l - iz) + c(l + iz)] = c + 1 

g{z) = Im[(l — iz) + c(l -t- iz)] = cz ~ z 

[A+(l + c)-i(a + 6)] =0,A = -(l + c)-i(a + ^) 
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F{z, y) = Re[(iy + a)(l - iz) + iyc(l + iz) + &(1 + iz)] = (1 - c) + a + 6 
G(z, y) = Im[(iy + a)(l - iz) + iyc{l + iz) + 6(1 + iz)] = y (c + 1) + z (6 - a) 
Obviously, /(z) = and .g(z) = have no common real root. 

By the condition (ii) in theorem 2, we have A = — (l + c)^^(a + &) < 0, that is (l + c)(5 + a) > 0. 
Finally, it is easy to see that F(z, y) = and G(z, y) — Q has no common real roots (z, y), z Q 
R,y e R~ {0}, if and only if (1 — c)(6 — a) < 0. Thus by theorem 2, we complete the proof. 
Corollary 2 The system 



x{t) + cx{t - 2t) + ax{t) + ax{t - r) = {1 + c ^ 0) 
is delay-independent stable if and only if 

a > and — 1 < c < 



(24) 



1 



:(l - z2) - z2 + 1 



Proof. Let 

/(z) = Re[(l-iz)2 + c(l + iz)2] 
g{z) = Im[(l - iz)2 + c(l + iz)2] = 2cz - 2z 
[A - (1 + c)-i(-a + -a)] = 0, A = -2a(l + c)-i 

F{z,y) — Re[(jy + a)(l — iz)^ + iyc(l + iz)^ + a(l + jz)(l — iz)] = a + 2yz — 2cyz + az'^ + a (l — z^) 
G(z, y) — Im[(iy + a)(l — iz)"^ + jyc(l + iz)^ + a(l + iz)(l - iz)] = y (l - z^) — 2az + cy (l - z^) 
Obviously, /(z) = and g{z) — have no common real root. 

By the condition (ii) in theorem 2, we have A = — 2a(l + c)^^ < that is a > 0, 1 + c > 0; or 
a < 0, 1 + c < 0. 

Finally, the resultant of F{z,y) and G{z,y) 
2z - 2cz 2a 



R{F,G){z) det 



(l-z2)(l + c) 



-2az 



= 2a(3c- l)z^ - 2a(l + c) 



It is easy to see that F{z, y) = and G(z, y) = has no common real roots (z, y) , z ^ R, y ^ 

R — {0}, if and only if a > and —1 < c < -. Thus by theorem 2, we complete the proof. 

Remark 8 Corollary 1 and Corollary 2 are consistent with the relevant results in I22| . 
respectively. They have also discussed delay-dependent stability of System and System 

Next, we will present some numerical examples using the main results of this paper. 

Example 1 Consider the following system 



X{t) - CX{t -t) = AX(t) + BX{t - t) 



where 



G 



0.1 




-3 -2 




1 






0.1 


1 


1 



and a is a nonzero constant. 

We now consider the maximal stability bound in terms of a. Let / be the 2x2 identity matrix, 
/(z) := Rc[dct[/(1 - iz) - G(l + iz)]] = 0. 81 - 1. 21z2 
g{z) := lm[det[/(l - iz) - G(l + iz)]] = -1. 98z 

det[A/ - (/ - C)-^{A + B)] = A^ + 3. 3333A + 2. 4691 + 1. 2346a - 1. 2346a2 

F(z,y) := Y{.c[Act[{iyI-A){l-iz)-iyC{l + iz)-B{l + iz)]] = y^ (l.21z2 - O.Sl) +6.0yz + a-2 

z"^ - + z'^a + z'^a^ + 2 

G(z,y) lTa[Aet[{iyI - A){l~iz) ~iyC{l + iz) - B{l+iz)]] = 1. OSy^z + y (2. 7 - 3. 3z2) -4z- 

2za2 

Obviously, /(z) and g{z) have no common real roots. 
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If ReA[(/ - C)~^{A + B)] < 0, by the Hurwitz criterion, we have -0. 99997 < a < 2.0. 
Finally, the resultant of F{z,y) and G{z,y) 

R{F, G) (z) (-5. 9049a + 5. 90490^ - 11.81) + (-57. 802 - 24. 592^2 _,_ gQa^ - 14. 256a3 _ 
11. 875a)z2 + (6. 336a-31. 68a''-81. 081a2-106. 01 + 3. l68a^)z'^ + {25. 483a-40. 6080^-86. 346+ 
17. 424a3 + 19. S6a^)z'^ + (13. 177a + 13. 177a2 - 26. 354)^8 

By a careful calculation, R{F,G){z) = has no real root when a G (—0.99997,1]. In fact, 
if a G (—0.99997,1], all coefficients of R{F,G){z) have the same sign. Thus, F{z,y) — and 
G{z, y) = have no common real roots. Hence, by Theorem 2, when a £ (—0.99997, 1], the system 
is delay-independent stable. 

Remark 9 The criteria in |13II19| do not work here. [21] concluded that the system is delay- 
independent stable when ]a] < 0.989. By using the criteria in this paper, we conclude that the 
system is delay-independent stable when a G (—0.99997, 1], which gives a less conservative bound 
of a. 

Example 2 Consider the following system 

Xit) - CX{t - t) = AX{t) + BX{t - t) 

where 



0.4 




-1 




1 


,B = a 


0.4 


-1 


1 



and a is a nonzero constant. 

We now consider the maximal stability bound in terms of a. Let / be the 2x2 identity matrix, 
/(z) := Re[det[/(1 - iz) - G{1 + iz)]] = 0. 84 - 0. 84z2 
g{z) Im[det[/(1 - iz) - C(l + iz)]] = -2. 32z 

det[A/ - (/ - C)-\A + B)]] = A2 + 2. 381A - 0. 95238Aa + 1. 1905 - 1. 1905a2 

F{z,y) := Re[det[{iyI-A){l-iz)-iyC{l + iz)-B{l + iz)]] = -0. 84y2 + 1 + 4zy + 0. 84z V - 
z^ + 1. 6yaz — a^ + a^z^ 

G{z, y) Im[dct[(i2;/ - A){1 - iz) - iyC{l + iz) - B{1 + iz)]] = 2y + 2. 32zy'^ - 2z - 2z'^y - 
0. 8ya + 0. Sz'^ya - 2a'^z 

Obviously, f(z) and g{z) have no common real roots. 

If ReA[(/ - C)~^{A + B)] < 0, by the Hurwitz criterion, we have -1.0 < a < 1.0. 
Finally, the resultant of F{z,y) and G{z,y) 

i?(F, G){y) := (-4.0 + S.Oa'^ - 4.0a®) + (2. 56a^ - 2. 56a6 - 16.0 + lQ.Qa'^)y^ + (7. 3856a4 + 
5. 12a2 - 23. 795)?/^ + (2. 4945a2 - 15. bQ)y^ - 3. 7978y® 

It is easy to see that R{F,G){y) — has no real root when a £ (—1.0,1.0). In fact, if 
a G ( — 1.0,1.0), all coefficients of R{F,G){z) have the same sign. Thus, F{z,y) = 0,G{z,y) = 
have no common real roots. Hence, when a G (—1.0, 1.0), the system is delay-independent stable. 

Remark 10 The criteria in [BUT^OU] can work here, and they can get a's bounds respectively 
as: Ja] < 0.2 [TS|; |a| < 0.2 [El; Ja] < 0.9165 EOI- The result of this paper is a G (-1.0,1.0), 
which obviously gives a less conservative bound of a. 

Example 3 Consider the following system 

X{t) - CX(t - r) = AX{t) + BX{t - r) 

where 



C = 



" 0.2 




"-20" 




" 0.5 " 






,B = 




0.2 


-1 


0.5 
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We will discuss whether the system is delay-independent stable. Let / be the 2x2 identity 
matrix. 

f{z) := Re[det[/(1 - iz) - C(l + iz)]] = 0. 64 - 1. Uz^ 

g{z) Im[det[/(1 - iz) - C(l + iz)]] = -1. 92z 

det[A/ - (/ C)-^{A + B)]] = A2 + 3. 75A + 2. 7344 = 

F{z, y) Y{.e[Aet[{iyI-A){l~iz)~iyC{l+iz)~B{l+iz)]] = -0. 64?/2+6.0z?/+l. 75+1. 44zV- 
1.75z2 

G{z,y) := lm\dci[{iyl - A){l-iz)-iyC{l + iz)- B{l + iz)]] = 2. 4y + 1. 92zy2 - 4. 5^ - 3. Gz^y 

Obviously, f{z) and g{z) have no common real roots. 

By the Hurwitz criterion, it is obvious Re\[I - (/ - C)^^{A + B)] < 0. 

The resultant of F{z, y) and G{z, y), R{F, G){y), is 

R{F, G){y) := -35. 684?/6 _ 122. 57^^^ - 3. 3974y8 _ 152. A6y^ - 62. 016 

R{F, G){y) — has no real roots, that is, F{z, y) — 0, G{z, y) = have no common real roots. 
Hence, the system is delay- independent stable. 

Remark 11 By the criteria in fQEOli the maximum allowable bound guaranteeing asymptotic 
stability of the system was obtained, they are Tmax = 0.1352 in |18| and Tmax = 0.7516 in |2U| . 
respectively. However, this paper concludes that the system is delay-independent stable. The 
reason for the different conclusions is that the criterion in this paper is necessary and sufficient, 
whereas the criteria in |18l I20j are only sufficient conditions. 

The above three examples were all discussed in Sharper stability bounds are obtained by 
using the algebraic criteria in Theorem 2. 

Example 4 Consider the system 

X{t) - CX{t -t) = AX{t) + BX{t - t) 



where 



0.5 








0.2 








0.3 







.B = 



-1 1 







We will discuss whether the system is delay-independent stable. Let / be the 3x3 identity 
matrix. 

f{z) Re[det[/(1 - iz) - G{1 + iz)]] = 0. 28 - 3. 6z^ 
g{z) := Im[det[/(1 - iz) - C(l + iz)]] = -1. 78z + 2. 34^3 
det[A/ - (/ - G)-\A + B)] =X^ + 16. 036A2 + 87. 857A + 160. 71 

F{z, y) := Rc[dct[(iyI-A){l-iz)-iyC{l+iz)-B{l+iz)]] = 45+72. 4zy+24. 75zV-l- 78y^z+ 
2. 342^2/3 _ 10. Az^y - 4. 49y2 - 35z'^ 

G{z,y) lm[det[{iyI-A){l-iz)-iyC{l + iz)-B{l + iz)]] = 19.89y'^z-52.6z'^y + 3.6y^z'^ - 
6. 87.2^^2 _ 28y3 + 5^3 _ 75^ 24. Qy 

It is easy to see that f{z) — 0, g{z) = have no common real roots. 

By the Hurwitz criterion, all roots of det[A/ — (/ — G)^^{A + B)] have negative real parts. 
Furthermore, by a careful calculation using Lemma 8, F{z, y) and G{z, y) have no common real 
roots. 

Hence, the system is delay-independent stable. 

Remark 12 Let C = in this example, then it becomes a retarded differential system. In this 
case, the conclusion in this paper is still correct and is consistent with the one in [71 130j . 

Finally, we present a more complex example on computing the maximal allowable delay bound. 
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Example 5 Consider the following system 



X{t) - CX{t - r) = AX{t) + BiX{t - t) + B2X{t ~ 2r) + B:iX{t - 3r) 



where 



C 



0.02 




-0.05 0.005 





0.01 





0.03 



0.5 








0.25 



.A 



0.005 


-1 



0.005 







0.25 





-0.5 



0.0375 


0.05 




.Bo = 





0.05 
0.05 
-2.5 



-2 -3 -5 
0.005 0.0025 





1 

-2 






_ -1 

0.075 

0.05 










-0.5 
0.125 ■ 



-1 



0.05 


-0.5 






0.0005 




Let / be the 4x4 identity matrix. 

f{z) := Re[det[/(1 - iz) - C(l + iz)]] = 0. 36383 - 5. 7048z2 + 1. 93162^ 
g{z) := Im[det[/(1 - iz) - C(l + iz)]] = -2. 4477z + 5. 5521z^ 

det[XI - (/ - C)^i (A + Bi + ^2 + S3)] = A-* + 3. 949A3 + 16. 661X^ + 19. 9A + 11. 646 
F{z,y) Re[det[{iyl - A){1 - iz)^ - iyC{l + iz){l - iz)'^ - Bi{l + iz){l - iz)'^ - ^2(1 + 
2(1 - iz) - Bsil + iz)^]] = 365. 27y2z2 + 4587. 5y'^z^ - 771. 21z^y - 2331. 4zJy 



911. 59?;V + 625. 81yS^ - 4. 821ziV 



35. A73y^z^ + 368. 62y^z* 



1.9316ziV 
26.68z"y3 _ 104. 21ziV 



2.008zi2_21.357z"y- 



-321. 27 z^"y^ + 1 189. 52^ - 2558. Ay^z^ + 64. 8A5zy + 2334. 2z^y - 6. 0618?/2 - 



211. 53^2-1765. 12^-2396. 7zV-555. 07zV+621. 06z^y-17. 459y^z+335. 39^3^3-1413. Gy^z'^ 



1747. 8zV 

G(z,y) : 



813. 91z«- 111.95^10 
= Im[det[(i?;/ - ^)(1 - 



.36383y4 + 4.2371 

j)3 - iyC{l + iz){l - - Bi{l + iz){l - iz)^ - ^2(1 + 



iz)^{l~iz)~ Bsil + iz) 



-5. 3583y'*z + 140. ly V + 1835. 7y^z'^ - 311. 72y V + 97. 231y 



4002. 2?;''z'' - 805. 43?/^z' 
682.09y'^z^ - 1171.7^^^3 



1542. 3z4?/-59. 037ziy-278. 51z'^y-163. Olz 



10, 



1457. 7z^y - 1190. Iz^y^ - 3875. Iz^y^ 



69. 443y2z + 0. 66541zi2y + 887. 73zV - 1- 6155zi2y3 + 21. 005z"?/^ 



?/ + 7.2403?/-44.703z- 
1062. 4zV - 2694. bz^y + 
164. 45ziV + 22. 407z" - 
1. 4368?;3 + 605. 71z3 - 1683. 4z-^ + 1381. 2z^ - 356. 38z^ 

Obviously, /(z) — 0,g{z) — have no common real roots. 

By Hurwitz Criterion, det[A/ — {I — C)~^ {A + Bi + B2 + B3)] has only negative-real-part roots. 

Further, by Lemma 8, we have the common real roots of F{z, y) = and G(z, y) = are {y = 
-0. 86798, z = 10. 823}, {y = 0. 58092, z = 10. 7}, {z = -1.6587,y = 1.3876},{y ^ 1.1338,z = 
-0. 42313}, {y = -0.58092, z = -10. 7}, {y = 0. 31594, z = -1.2099},{z = -10.823,?/ = 
0. 86798}, {y = -1. 1338, z = 0. 42313}, {y = -3. 3266, z = 0. 24369}, {y = 3. 3266, z = -0. 24369}, 
{y = -1. 3876, z = 1. 6587}, {y = -0. 31594, z = 1. 2099}. 

Therefore, the system is not delay- independent stable by Theorem 2. The maximal delay 
bound guaranteeing the stability is T = 0.14371 (by taking all common real roots of F{z,y) = 
and G{z,y) = into ^). 

Remark 13 Set C = in Example 5, then the system degenerates to a retard differential 
system. The maximal delay bound was obtained in as T = 0.4777. The same conclusion can 
be obtained by using the results in this paper. 
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From the above corollaries and examples, it is obvious that our criteria work well on judging the 
delay-independent stability of neutral differential systems, and its simplicity, accuracy, convenience, 
and wide applicability greatly facilitate the engineering practice. 

6 Conclusion 

This paper establishes some algebraic criteria for determining the delay-independent stability of a 
class of neutral differential systems x{t) — X^feLi Bkx{t — kr) = Aox{t) + ^kx{t ~ kr), and 

presents a method for determining the maximal delay bound guaranteeing stability if the systems 
are not delay-independent stable. To the best of our knowledge, the algebraic criteria has some 
noteworthy characters comparing with the related literature. 

• Generality 

The system x{t) — B]^x{t — fcr) = AQx{t) -I- X^fcLi — fcr) is more general, and covers 

various forms of systems studied in|Slia[niiniIISlIIlIIllSIIE21ESlEll2Sll2niEni- The criteria 
in this paper can be applied to all these systems, and the relevant results are consistent with or 
better than that of those literature, which have been illustrated by corollaries and examples. 

• Necessary-Sufficient Conditions 

The criteria for delay-independent stability are necessary and sufficient, and the broader con- 
ditions can be got comparing with that of some literature (see example 1-3). 

• Practicality 

The conditions in the criteria are all algebraic conditions. Condition (ii) can be checked by 
Hurwitz Criterion Determining real roots of polynomials in conditions (i) and (iii) can be 
carried out by the "Complete Discrimination System for Polynomials" mentioned before. The sign 
list of the discrimination sequence of polynomials with symbolic coefficients can be obtained easily 
by computer |27[ I28[l^ . Therefore, "on-line" operation can be realized. 

References 

[1] Avellar, C. E. and Hale, J. K., On the zeros of exponential polynomials, J. Math. Anal. AppL, 
1980, 73: 434-452. 

[2] Bellman, R. and Cooke, K. L., Differential-Difference Equations, Academic Press, New York, 
1963. 

[3] Boukas, E. K. and Liu, Z. K., Deterministic and Stochastic Time Delay System, Birkhauser, 
Boston, Basel, Berlin, 2002. 

[4] Brumley, W. E., On the asymptotic behavior of solutions of differential-difference equations of 
neutral type, J. Diff. Eqs., 1970, 7: 175-188. 

[5] Chen, J. and Latchman, H. A., Frequency sweeping tests for stability independent of delay, 
IEEE Trans. Automat. Contr., 1995, 40: 1640-1645. 

[6] Gantmacher, F., The Theory of Matrices, Chelsea, New York, 1959. 



14 



[7] Gu, N., Tan, M. and Yu, W. S., An algebra test for unconditional stability of linear delay 
systems, Proc. of the 40th IEEE Conference on Decision and Control (CDC 2001), Orlando, 
Florida, USA, 2001, 4746-4747. 

[8] Hale, J. K., Theory of Functional Differential Equations, Springer- Verlag, New York, 1977. 

[9] Hale, J. K., Infante, E. F. and Tsen, F. -S. P., Stability in linear delay equations, J. Math. 
Anal. AppL, 1985, 105: 533-555. 

[10] Hale, J. K. and Verduyn Lunel, S. M., Introduction to Functional Differential Equations, 
Springer- Verlag, New York, 1993. 

[11] Han, Q. -L., On delay-dependent stability for neutral delay-differential systems. Int. J. Appl. 
Math. Comput. Sci., 2001, 11: 965-976. 

[12] Hertz, D., Jury, E. J. and Zeheb, E., Stability independent and dependent of delay for delay 
differential systems, J. Franklin Inst, 1984, 318: 143-150. 

[13] Hu, G. Di and Hu, G. Da, Some simple stability criteria of neutral delay-differential systems, 
Appl. Math. Comput, 1996, 80: 257-271. 

[14] Kamen, E. W., Lineal systems with commensurate time-delays; Stability and stabilization 
independent of delay, IEEE Trans. Automat Contr., 1982, 27: 367-375. 

[15] Kharitonov, V. L., Robust stability analysis of time delay systems: A survey. Annual Reviews 
in Control, 1999, 23: 185-196. 

[16] Kharitonov, V. L., Torres-Muiioz, J. A. and Ortiz-Moctezuma, M. B., Polytopic families of 
quasi-polynomials: vertex-type stability conditions, IEEE Trans. Circuits Syst., Part I, 2003, 
CAS-50: 1413-1420. 

[17] Kharitonov, V. L. and Zhabko, A. P., Robust stability of time-delay systems, IEEE Trans. 
Automat Contr., 1994, 39: 2388-2397. 

[18] Khusainov, D. Ya. and Yun'kova, E. V., Investigation of the stability of linear systems of 
neutral type by the Lyapunov function method, Diff. Uravn., 1988, 24: 613-621. 

[19] Li, M. L., Stability of linear neutral delay-differential systems. Bull. Austral. Math. Soc, 
1988, 38: 339-344. 

[20] Park, J. -H. and Won, S., Stability analysis for neutral delay- differential systems, J. Franklin 
Inst, 2000, 337: 1-9. 

[21] Qin, Y. X., Liu, Y. Q., Wang, L., and Zheng, Z. X., Stability of Dynamical Systems with 
Delays (Second Edition), Science Press, Beijing, 1989. 

[22] Ren, H. S. and Li, H. Y., Explicit asymptotic stability criteria for neutral differential equations 
with two delays, Appl. Math. E-Notes, 2002, 2: 1-9. 

[23] Ren, H. S. and Zheng, Z. X., Algebraic criterion for asymptotic stability of neutral equations 
of the form x{t) + cx{t -t) + ax{t) + bx{t - r) = 0, Acta. Math. Sinica, 1999, 42(6): 1077-1088. 

[24] Verriest, E. L and Niculescu, S. I., Delay-independent stability of linear neutral systems: A 
Riccati equation approach, in: Stability and Control of Time-Delay Systems (eds, L. Dugard 
and E. L Verriest), Springer- Verlag, London, 1997, 92-100. 



15 



[25] Wang, Z. D., Lam, J. and Burnham, K. J., Stability analysis and observer design for neutral 
delay systems, IEEE Trans. Automat. Contr., 2002, AC-47: 478-483. 

[26] Wu, J. Y., The position of boundary surface for a sort of trancendental equation with three 
parameters. Acta. Math. Sinica, 1994, 37: 301-308. 

[27] Yang, L., Recent advances on determining the number of real roots of parametric polynomials, 
J. Symbolic Computation, 1999, 28: 225-242. 

[28] Yang, L., Hou, X. R. and Zeng, Z. B., A complete discrimination system for polynomials. 
Science in China, 1996, E-39: 628-646. 

[29] Yang, L., Zhang, J. and Hou, X., Nonliear Algebraic Equations and Machine Proving, Shanghai 
Science and Education Press, Shanghai, 1996. 

[30] Yu, W. S. and Wang, L., Algebraic criteria for dclay-indcpcndcnt stability of differential system 
x{t) = Aox{t) + Y.^^^ AKx{t - kr), Chinese Control Conference (CCC2000), Hong Kong, 2000, 
435-439. 



16 



